Exact solution of the Faddeev-Volkov model 
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The Faddeev-Volkov model is an Ising-type lattice model with positive Boltzmann weights where 
the spin variables take continuous values on the real line. It serves as a lattice analog of the sinh- 
Gordon and Liouville models and intimately connected with the modular double of the quantum 
group U q {sl2). The free energy of the model is exactly calculated in the thermodynamic limit. In 
the quasi-classical limit c — > +00 the model describes quantum fluctuations of discrete conformal 
transformations connected with the Thurston's discrete analogue of the Riemann mappings theo- 
rem. In the strongly-coupled limit c — > 1 the model turns into a discrete version of the T> = 2 
Zamolodchikov's "fishing-net" model. 
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Faddeev and Volkov [l|, 0, El obtained a very interest- 
ing solution of the star-triangle relation. There is only 
a few two-dimensional solvable lattice models [4j where 
the Yang-Baxter equation takes its distinguished "star- 
triangular" form. These models include the Ashkin- 
Teller [j| , Kashiwara-Miwa Q and chiral Potts models 
(the latter also contains the Ising, self-dual Potts @] and 
Fateev-Zamolodchikov [j| models as particular cases). 
All above models are also distinguished by a positivity 
of Boltzmann weights — the property that is naturally 
expected for physical applications, but rarely realized for 
generic solutions of the Yang-Baxter equation. Recently, 
we have observed [lfj that the Faddeev-Volkov model 
[E 0, 0] also possesses the positivity property. Apart 
from being an interesting model of statistical mechanics 
and quantum field theory in its own rights (it serves as 
a lattice version of the sinh-Gordon and Liouville mod- 
els [ll, 12 [) this model has remarkable connections with 
discrete geometry. As shown in (lCj| it describes quan- 
tum fluctuations of circle patterns [13J and associated 
discrete conformal transformations connected with with 
Thurston's discrete analogue of the Riemann mapping 
theorem 
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FIG. 1: The square lattice (solid lines) and its medial 
"rapidity" lattice (dashed lines). 



Consider the square lattice, shown in FigfTJ Each site i 
of a lattice is assigned with a spin variable ai £ R, taking 
continuous real values. Two spins a and b interact only if 
they are adjacent (connected with an edge of the lattice). 
Typical horizontal and vertical edges are shown in FigJ^l 
The corresponding Boltzmann weights are parametrized 
through the so-called "rapidity variables" p and q associ- 
ated with the oriented dashed lines in FigJT] In our case 
the weights depends only on the spin and rapidity differ- 
ences a — b and p — q, where a and b are the spins at the 



ends of the edge. We will denote them as W p 



b) 



and W p - q (a — b) for the horizontal and vertical edges, 
respectively. The partition function is defined as 



Z = 



Y\_W p - q (<Ti-cTj) Y\_W p - q (a k - ai) Y[ d(J i 



(1) 

where the first product is over all horizontal edges 
and the second is over all vertical edges (k,l). The in- 
tegral is taken over all the interior spins; the bound- 
ary spins are kept fixed. The weights W p - q (a — b) and 
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FIG. 2: Two types of Boltzmann weights. 

Wp-q(a_— b) are related with non-compact representa- 
tions 15] of the modular double [3] of the quantum 
group U q (sl2)®Uq(sl2), where q = e l7rb and q = e~™/ b . 
The modular parameter b is connected to the Liouville 
central charge Ci = l + 6(b + b -1 ) 2 . It is convenient to 
define 



exp 



i7T(b-b- 1 )' 

2(b + b"l), 



(2) 
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The main physical regimes of the model 

(i) b > 0, and (ii) |b| = 1, Im(b 2 ) > 0. (3) 

correspond to real values of rj. For the regime (i) it is 
enough to consider the range < b < 1 (due to the 
symmetry b <-* b _1 ). 

Introduce the non-compact quantum dilogarithm 



ip{z) 



exp 



-2\zw 



dw\ 



4 7 R+i0 sinh(wb)sinh(w/b) w J 
-qe 2 " b ; q 2 )oo 



(-qe 



';q 2 )c 



(4) 



(1 — q 'x). The product represen- 

k— 

tation above is valid for Im(b 2 ) > 0. 

The Boltzmann weights Wg (s) and Wg (s) are defined 

as 

^)4^t^7T' We(s) = W«Ms), 
F e tp(s - \r]9 I tt) 

(5) 

where 9 and s stand for the rapidity and spin differences, 
respectively. The normalization factor Fg has the form 



F u 



where 



$(z) 



ex P( - 

K+iO 

(qV^q 4 )^ (-q e -/^;q 2 ) 



(6) 



dw\ 



sinh(uib) sinh(ui/b) cosh(2ui77) w ' 

(7) 
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with q defined in (JU). The main properties of (f{z) and 
<f>(z) are discussed in (To| . 

When the parameter b belongs to either of the regimes 
(J3|), the Boltzmann weights Wg (s) and Wg{s) are real 
and positive for < 9 < 7r and real s. They are 
even functions of the variable s and decay exponentially, 
Wg(s) ~ F 9 ~ 1 e~ 2r ' e l s l, when s -> ±oo. The weight 
Wg(s), considered as a function of 9 at fixed real s, is 
analytic and non-zero in the strip < Re 9 < tt . 

The weights Wg and Wg satisfy several functional re- 
lations: the duality relation 



Wg (s) 



the inversion relations, 



dx e 



2tt\xs 



Wg (x) , 



(8) 



lim / dcWit+ e (a-c)W- lt +e(c-b) = S(a-b) 

E--0+ .' 



We (a - (a - 6) = 1 , (9) 



where t is real, and the star-triangle relation, Fig(3j 
da W q ^ r (a — a) W p ^ r (c — a) W p - q (a — b) 

= W p . q (c - a)W p - r (a - b) W q . r (c-b). (10) 





FIG. 3: Star-triangle relation. 



We used the inversion relation method jl7, 18, to 
exactly calculate the partition function ([T]) in the thermo- 
dynamic limit. The result is included in the normaliza- 
tion of the of Boltzmann weights, so that the free energy 
per edge, 



(3 f edge = - lim TV^log^O 



(11) 



vanishes when the number of edges, N, tends to infin- 
ity (provided the boundary spins are kept finite). The 
weights (O attain their maximal values at s = 0. In the 
quasi-classical regime (i), < b < 1, the value of the 
constant Wn. (0) slowly interpolates between the values 



= e-- 



b=l 



V2, (12) 



where G = 0.915965 ... is the Catalan's constant. The 
s-dependence of the two-spin interaction energy E(s) = 
— \ogW^(s) is quadratic, E(s) — E(0) ~ s 2 , for small s 
and gradually becomes linear, E(s) ~ 7rr/|s|, for large s. 

The parameter b 2 plays the role of the Planck constant. 
The quasi-classical limit b — > (where the model reveals 
remarkable connections with the discrete geometry) was 
considered in [Toj . It is worth mentioning some details 
here. When b — > the weight function Wg(s) acquires a 
very narrow bell-shaped form and rapidly decays outside 
the small interval \s\ < b/9. The partition function ((T|) 
can be then calculated with saddle point method. It is 
convenient to define 9 = p — q and use rescaled spin 
variables p — {pi, p%, . . .}, given by pi — 2irb<Ji. The 
leading asymptotics of (J5]) reads 



T4Vp/27rb) = e^l^ b2 + ( b °) 







where the function 
rP 



A{9\p) = - /'iGgMOde, fe(0 
1 Jo 



1 + e^+ ie 



(13) 



(14) 
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is simply related to the Euler dilogarithm. The station- 
ary point of the integral (p} is determined by the classical 
equation of motion 

fe{pw-p)U-e(pN-p)fe(pE-p)U-e(ps-p) = 1, (15) 

connecting five neighboring spins pw, Pn, Pe, PS and p, 
arranged as in FigU This non-linear difference equation 
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FIG. 4: Arrangement of the p- and i/j-variables around 
a four-edge star. 

is a particular (square lattice) case of the "cross-ratio 
equations" [UEH- The latter can be viewed as discrete 
analogues of the Cauchy-Riemann conditions for analytic 
functions. In particular, Eq . (| 1 5|) defines discrete confor- 
mal transformations [2y| of the square lattice [2(| l22l |. 
where the local dilatation factors r$ = exp(p.j) are de- 
termined by the rescaled spin variables p.;, solving (|I5p . 
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sites. 

The cross ratio equations (TT5]) are closely related to the 
Hirota difference equation [231 ] . Place a purely imaginary 
variable tpi, < Im^ < 2tt, on every site i of the dual 
lattice. These sites are shown by black dots located in 
Figs. 0] and [5l Connect these new variables to the existing 
variables pi by the following rules. Let iptf and ips be the 
variables located above and below of a horizontal edge, 
as in Fig. [5) Then, we require that 



qPW _i_ qPe+W 
qPe _|_ qPw+W 

Similarly for a vertical edge, 

qPn _l qPs+i8" 



(16) 



^IpE-lpW 



gPS -|_ qPn+iS' 



q4>w _i_ q4>e+>8 

(17) 



where 9* = tt — 9. The consistency of these definitions 
across the lattice is provided by (| 1 5|) . Note that the sec- 
ond form of (|17p is identical to ptj|) upon interchanging 
all tp- and p- variables. Therefore it is natural to associate 
this universal equation with every face of the "double" 
lattice consisting of all white and black sites (its faces 
are shown in Fig. [5|- This is the famous Hirota equation 
introduced in [23| . 

Thus, in the quasi-classical limit b — ► the stationary 
points of the integral (pj defines discrete analogs of con- 
formal transformations. At finite values of b the model 
describes quantum fluctuations of these transformations. 
Given that the spins pt define the local dilatation factors 
the Faddeev-Volkov model describes a quantum discrete 
dilaton. The continuous quantum field theory with the 
conformal symmetry (24| has remarkable applications in 
physics and mathematics. It would be interesting to un- 
derstand which aspects of the continuous conformal field 
theory can be transferred to the discrete case. 

For the "strongly coupled" regime (ii), when the pa- 
rameter b is on the unit circle |b| = 1, the function Wtl (s) 
still has its absolute maximum at s = but starts to 
develop additional side maxima at integer points when 
arg(b) approaches the value tt/2. The first pair of such 
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FIG. 6: The function W 7T / 2 {s) for T) = 0.30 and 0.06. 

side maxima at s = ±1 appears when rj ~ cos(27r/5) ~ 
0.309, see Fig. [6l For rj = the function W*(s) turns 
into a superposition of (5-function like peaks at ,s G Z. In 
particular, the height of the central maximum diverges 
as 

(o) ~ (2^r 1/2 r(|)/r(|) + 0(^/2), n - o . 

(18) 

A more careful limiting procedure capturing a fine struc- 
ture of these sharp peaks requires a redefinition of the 
spin variables. Using the asymptotics 

(p(n + xr)) ~ e~ i7r / 12 (47T77) i3 T( 1 - n + ix ) /r( h^s. ) ; 

$(2x7?) ~ e-' ,7T / 24 (8T:r 1 y x r(^)/r(^), (19) 
where rj — > 0, n € Z and |x| -C f] 1 ' it is easy to see that 
We{n + xrf) ~ 7]- 6 ^V e {n, x), rj -> (20) 
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where the function 

Vg(n,x) 



1+8 /n 1 — n — 6 /■k+\x \^Y'^ l — n — 6/7T — \x ^ 



(2tt)^ r (iz|ZzL)r( 1 -"+^ 7r+i:E ; 



1 — 7i-\-9/tt — \X\ 

(21) 

is real and positive for n € Z and x £ R. This func- 
tion defines Boltzmann weights for a new model where 



J 



each lattice site i is assigned with a pair of fluctuating 
variables (rii,Xi), where ii; G Z take integer values and 
Xi £ K take continuous values on the real line. Its parti- 
tion function is defined by ([T]) where Wg(o~i — o~j) replaced 
by Ve(n,i~ rij , Xi — Xj ) (and similarly for W ) and every in- 
tegral over Oi is replaced by a sum over rii and an integral 
over Xi, namely J doi — * S n . eZ / dxj. 



The star-triangle relation for Vg (n, a;) simply follows from (I10|) , 

/ dxoVfiiCni - n ;xx - x )Vg 2 (ri2 - no; £2 - x a )Ve 3 (n3 - n ;x 3 - x a ) 

1 ^rj.JVL 



(22) 



= K--Si («2 -n 3 ;x2- aj 3 )K-e a («i - n 3 ; Xi - x 3 )K--e 3 ("i ~ n 2; xi - x 2 ) 
where 6>i + 6*2 + 8 3 = 2tt. Similarly, the inversion relations © imply 



lim V] / dx 2 K-+e(ni - n 2 ;xi - x 2 )K--e(^2 - n 3 ;x 2 - x 3 ) = S ni , n3 5(xi - x 3 ) , Vg(n,x)V-g(n,x) = 1 (23) 
- >0+ n a6 z- /R 



The Boltzmann weight Vg(n; x) is an even function of n 
and x; it satisfies the following initial conditions 



V n (n;x) = 6 n , 5(x) , V (n;x) = l. 
Note that its asymptotics when n,x — > ±00, 

•2^r(i±§^) . 2 , 2l _ e/7r 



Ve(n;x) ~ (-) 



r(i^) 1 

coincides with the normalized Boltzmann weight of the 
2> = 2 Zamolodchikov's "fishing-net" model [25| (except 
that the variable n in our case is discrete). It follows 
from (fTTj) that with this normalization the specific free 
energy of the "fishing-net" model vanishes in the ther- 
modynamic limit in complete agreement with 25]. 



(24) 



(35) 
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